Magnetoelectric effect in mixed valency oxides mediated by charge carriers 
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We show that the presence of free carriers in a substance can generate the multiferroic behavior. 
Namely, if the substance has mixed-valence ions, which can supply free carriers and have electric 
dipole and spin moments, all three types of long-range order (ferromagnetic, ferroelectric and mag- 
netoelectric (ME)) can occur at low temperature. The physical origin of the effect is that charge 
carriers can mediate the multiferroic behavior via spin - spin (RKKY), dipole-dipole and dipole 
- spin interactions. Our estimate of the interaction magnitude shows that there exist an optimal 
carrier concentration, at which the strength of ME interaction is maximal and comparable to that of 
spin-spin RKKY interaction. This permits to conclude that in substances, where RKKY interaction 
between local spins is not small, a substantial value of free carriers mediated ME interaction can 
occur. Our analysis shows that disorder in the above substances does not suppress multiferroic 
effects. 



PACS numbers: 77.80.-e, 77.84.Bw, 75.50.Dd, 77.80.Bh 

Multiferroics are substances which can simultaneously 
exhibit spontaneous polarization, magnetization and/or 
strain The ferroics where magnetic and electric or- 
ders coexist have recently become the focus of intensive 
research [1, H, H|, magnetism and ferroelectricity being 
involved with local spins and off-center structural distor- 
tions (electric dipoles) respectively. The physical nature 
of the coexistence of these two seemingly unrelated phe- 
nomena as well as of the recently observed Q strong in- 
teraction between polarization and magnetization (mag- 
netoelectric (ME) coupling) is not completely understood 
up to now, because trivial spin-orbit interaction can be 
the source of a very small ME effect only. The clarifi- 
cation of the physical nature of above phenomena would 
be of great importance both for fundamental research 
and for creation of new generations of electronic devices. 
In particular, the ME effect in multiferroics opens the 
way to write an information with the help of the electric 
field and to read it with the help of the magnetic field 
making the memory devices more stable and reliable Q. 
Among different multiferroics with high ME effect the 
oxides seem to be the largest group. Namely, perovskite 
rare-earth manganites RMnC>3 (R = Eu, Gd, Tb and 
Dy) HH, perovskite-typeBiFeOa, BiMn0 3 (see Ref. 
and references therein), materials with general formulae 
RMn 2 5 (R = Tb, Y, Eu) @, Q and other oxides be- 
long to this group It is easy to see that one of these 
oxides cations belongs to the ions of 3d-group known to 
have several valency states. Such ions are Mn 2+ , Mn 3+ , 
Mn 4+ , Fe 2+ , Fe 3+ , Cr 3+ , Cr 5+ etc. The presence of these 



ions have been observed by ESR and optical methods [§] 
so that mixed valency is expected to occur in the above 
materials. Moreover, due to chemical and structural com- 
plexity of RMn 2 5 , Mn 4+ ions are octahedrally coordi- 
nated by oxygen, whereas Mn 3+ are at the base center of 
a square pyramid 0] so that the general formula of this 
oxide can be rewritten as R(Mn 4+ Mn 3+ )Os with equal 
number of Mn 4+ and Mn 3+ . So, the cation valency de- 
pends on the local symmetry, defined by the number of 
surrounding oxygen ions. However it is known that the 
oxides used to have a lot of oxygen vacancies (and possi- 
bly other imperfections) randomly distributed in the host 
lattice. This can lead to a random distribution of cations 
with different valency, some of them can be magnetic or 
nonmagnetic (if their ground state spin is zero) . On the 
other hand, oxygen vacancies can generate free carriers so 
that the aforementioned oxides can acquire nonzero con- 
ductivity equal to that of semiconductors. We emphasize 
here, that this finite conductivity does not mean complete 
destruction of possible ferroelectric order, rather, many 
interesting effects appear, see for details. These ef- 
fects, however, are not directly related to the ME effect 
mechanism, discussed here. 

As we have mentioned above, the cation valency deter- 
mines its position either in the center of high symmetry 
surrounding or in the off-center position. This is because 
the position is defined by delicate balance of repulsive 
and attractive interactions between the ion and its sur- 
roundings, so that smaller ions would be expected to have 
the energy minimum in the off-center position (see Ref. 
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[3| and references therein). 

In Ref. [H, the data related to off-center positions of 
cations in spinel structure were reported. Note, that off- 
centrality of an ion can__be considered as a result of a 
pseudo-Jan Teller effect 
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It had been shown [l4( that in diluted magnetic semi- 
conductors (DMS) (for instance Gai-^Mn^As), also be- 
longing to above class of mixed valency substances, fer- 
romagnetism appears due to Ruderman-Kittel-Kasuya- 
Yosida (RKKY) indirect coupling between Mn d-shell lo- 
calized magnetic moments mediated by the induced spin- 
polarization in a free-hole itinerant carrier system. The 
experimental investigation of (Ga,Mn)As DNS [l5[ shows 
that the temperature of ferromagnetic phase transition is 
T c =110 K or even T c =173 K [la ], which is quite a large 
value allowing for the randomness of the Mn ions distri- 
bution. Since it is not excluded, that the ion (for example 
off-center Mn) can have both spin and electric dipole mo- 
ment, the above substances could be the "candidates" 
for realization of all three types of long-range order - 
ferromagnetic (due to RKKY interaction), ferroelectric 
(due to carrier mediated dipole-dipole interaction, see, 
e.g., Ref. [It]]) and ME order, also carrier-mediated. This 
means that the interaction via charge carriers can make 
the above substance a multifcrroic. Note, that phases 
with ME order belong to secondary ferroics with order 
parameter involving both magnetization and polarization 
[111 ] , their domain state can be switched by simultaneous 
application of magnetic and electric fields [12]. 

If local spins and electric dipoles (off-center ions) are 
randomly distributed over the host lattice sites, the over- 
all thermodynamic behavior depends on the history of 
the sample i.e. heating and/or cooling. These effects, 
characteristic for disordered systems [l8|, have been ob- 
served in ferroics with a high ME effect [2j. 

In the present paper we show that the microscopic na- 
ture of the magnetoelectric (ME) effect might be the in- 
teraction between local spins and dipoles via free charge 
carriers in a mixed valence oxide. In other words, if the 
substance under consideration has free carriers, they can 
mediate the ME effect. In the mean field approximation, 
we calculated for the first time the magnitude of carrier 
mediated ME interaction and compared it to both con- 
ventional RKKY interaction and dipole - dipole one. We 
obtain that the ME interaction oscillates with the ratio 
kpd (kp is a Fermi wave vector, d is off-center ion dis- 
placement) and find the optimal value of kp (and hence 
the free electrons concentration n e ) when the ME inter- 
action magnitude is the same as that of RKKY interac- 
tion. The dipole - dipole interaction magnitude is always 
smaller then both RKKY and ME. This permits to hope 
that in the substances, where the RKKY interaction be- 
tween local spins is not small, one can also achieve sub- 
stantial value of a ME interaction along with multiferroic 
state. 

The global symmetry of the linear ME effect implies 



the simultaneous breaking of time-reversal and inver- 
sion (say "space-reversal") symmetries of a substance. To 
break both above symmetries simultaneously, one needs 
to have in a system both spin moments S (axial vectors, 
changing their directions at time reversal) and dipole mo- 
ments d (in other words, two level system (TLS), see, 



e -g- ill) i changing their directions at spatial coordi- 



nates direction reversal. For a microscopic description 
this means that a system Hamiltonian should incorporate 
the "usual" terms proportional to J^ij Jss{rij)SiSj and 

y)jj Jddifi^didj (responsible for magnetic and electric 
long-range order respectively) as well as coupling terms 
proportional to - Jsd{'r f ij)Sidj. In our model, the lat- 
ter term is responsible for the ME effect in a solid. 

To account for above effect microscopically, we 
should find explicitly the potentials of interaction Jab 
(A,B=S,d). The ubiquitous method of calculation of the 
above interaction potentials Jab is a perturbation the- 
ory. In this approach, we have two " zeroth" Hamiltonians 
- free electrons, TLS'ses, or spins. The interaction term 
might be dipole-dipole (dd), spin-spin (SS) or dipole-spin 
(Sd) interaction. If this term is considered to be small, 
we can apply perturbation theory and get the interaction 
potential, as usually, in its second order. The details of 
calculations of dd and SS interactions can be found in 
Refs. 17, [3] and [2^, 21] respectively. Our aim is to 



calculate Sd interaction. This will be done similarly to 
above interactions, using the second order correction to 

(2) 

the ground state energy 

(A,B=S,d). This expression 

reads [H 

?(2) -J2 J ^ SiSk ' 
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where for finite temperatures 
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Here C/(ki — k 2 ) is the corresponding interaction ampli- 
tude, /i is a chemical potential. Actually, the specific 
form of J(r) (i.e. dd, SS or Sd) depends on the form of 
C/(ki— k 2 ). Here £(k) is a dispersion law for free electrons 



2 1,2 



2m* 



(3) 



and ni^ are Fermi-Dirac factors. 

Although the finite temperature effects might be im- 
portant, it can be shown that all important physics comes 
from the zero temperature case. This means that we may 
safely use the zero temperature approximation, which is 
actually used, for example, in RKKY interaction deriva- 
tion [211 ]. At zero temperature the occupation num- 
bers are simply unit step functions 0(£(k) — e/), where 
fx(T = 0) = 6f is Fermi energy. 
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The coupling function U(ki — k 2 ) is expressed via the 
expansion over spherical harmonics of the angle #ki,k 2 
between vectors \ti and k 2 , see Ref. [l7| for details. It 
can be shown that for our purposes it is sufficient to 
consider the function U to be isotropic, i.e. to use the 
zeroth term of this expansion 



I7(ki - k 2 ) = a 0(ki - k 2 ) 



(4) 



and the function <f>(ls.i — k 2 ) actually determines the kind 
of interaction (SS,dd or Sd). 

The spin-spin term J55 is indeed a well-known RKKY 
interaction, in this case </>ss(ki — k 2 ) = i and 



J(r) 



mV 2 alk i f sinx- 



-,x 



2k f r, 



(5) 



where V is a crystal volume. 

For calculation of dipole - dipolc coupling the in 
teraction function should be chosen in the form 17 
0dd(ki — k 2 ) = i sin [^d(ki — k 2 )] . Substitution of this 
function into Eq. ([2]) at zero temperature after a straight- 
forward but lengthy calculations yields 



J(r) 

F(y) 
a± 



mV 2 a?.kt 

° rT l[2F{x)-F^ + )-F^.)] 



8tt 3 ^ 2 
sin y — y cos y 
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x = 2kfr, ip± = a±x, 



v / l + A 2 ±2Acosa, A=-,d=|d|, (6) 



a is the angle characterizing the direction of vector d (for 
example, it can be thought as the angle between d and 
z axis. It is well-known (see above) that the off-center 
ions can have several orientations in a host lattice. The 
different orientations mean different values of angle a. 
For example, for a two position orientable dipole a = 
0, 7r. This fact will be used below. 

The Sd interaction potential will also depend on the 
angle a, its function |C/(ki — k 2 )| 2 = ctQ^gg^dd + c.c.) 
(c.c. stands for complex conjugated). This yields 



Jsd(r) 



47T 3 fr 2 

r)± = b±x, b± = J 1 



[F( V -) - F( v+ )], 



1 



A 2 ± A cos < 



(7) 



The expression ([7]) gives the desired potential of spin- 
dipole interaction mediated by charge carriers. 

To estimate the magnitude of the ME interaction, we 
calculate Weiss mean fields, related to SS, Sd and dd in- 
teractions. Having the interaction Jab(i") (A, B — S,d), 
the Weiss field is defined as 



W A B = yj J AB (r)d 3 r, 



(8) 
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FIG. 1: The magnitudes of SS, Sd and dd interactions (Weiss 
fields) as functions of parameter kpd at different directions of 
the off-center shift, defined by angle a. 



where V is the crystal volume. Since all Jab(i") except 
RKKY ([5]) depend on d - the ion off-central shift, the 
resulting integrals will be oscillating functions of this pa- 
rameter (to be more specific, the parameter is kpd). This 
means that depending on kpd value, the corresponding 
Weiss fields may be positive or negative, defining the mu- 
tual direction of spins and dipoles in ferroelectric and 
magnetoelectric phases. 

We recollect here, that in a mean field approximation 
the temperatures of the transition to ferromagnetic, fer- 
roelectric and magnetoelectric phases are related to cor- 
responding Weiss fields in an obvious fashion T c ab = 
y^nAnsl Wab\, where ua.b is a concentration (number 
of particles per unit volume) of local spins or electric 
dipoles. 

The expressions (|6]) and ((7]) show that the correspond- 
ing integrals cannot be done analytically. Substitution of 
the RKKY potential (0) into Eq. (HJ) yields 



W ss = \w , Wo 



m*Val 
n 2 fi 2 



(9) 



The estimate for Wq here gives Wq ~ 500K [17j. In 
that case Wss ~ 125K. The results of numerical calcu- 
lation of Wsd (ME interaction) and Wdd (dipole-dipolc 
interaction) are reported in Fig.l at different angles a. 
It is seen that while the magnitude of the magnetic in- 
teraction Wss is constant, the magnitudes of dd and Sd 
interactions oscillate with the parameter kpd. Moreover 
|VF d d| < |Wsd| < W ss and at a = and kpd » 1.4 
the value of |Wsd| is maximal, it is almost equal to 
Wss — 0.25VFo and is much more then |Wdd|< Actually, 
the ME effect can be realized in broader interval (see 
MD curve for a = on Fig.l) |FF S d| > 0.5|W SS | so that 
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0.5 < kpd < 1.4. This means that for a substance, hav- 
ing oxygen vacancies (which is a source of free carriers) 
and ions with off-center shift d there is strong possibility 
for substantial ME coupling via free carriers. 

We estimate carriers concentration in the above 
kp interval with the help of expression n e = 
Jn^ny eicpj—Ef/T), valid for low temperatures T < 100 
K[22|. where no ~ (10 22 — 10 23 )cm -3 is high temperature 
value of carrier concentration, ny rs 10 17 cm~ 3 is oxygen 
vacancies concentration, Ej — h 2 k'j/(2m*) is Fermi en- 
ergy and to* is carrier effective mass. For d — 0.5 ■ 10~ 8 
cm, m* = 10to € we obtain n e = (10 16 — 10 17 ) cm" 3 , 
which is a typical value for semiconductors. 

If the above substances have other (then oxygen va- 
cancies) sorts of defects, which can supply the free charge 
carriers, the above interactions are stronger, promoting 
the multiferroic behavior at low temperatures (T < 100 
K). Note, that many single-phase multiferroics exhibit 
multiferroicity only at low temperatures [2j. Our analy- 
sis shows (see Refs 23, [24j, |25| ) that the randomness in 
the above ion positions as well as presence of other types 
of unavoidable technological defects in a substance plays 
a dual role. On one hand, it tries to destroy all three 
kinds of above long-range order, the most sensitive fer- 
roelectric one, which is due to dipole-dipole interaction. 
On the other hand, this disorder might supply extra free 
electrons or holes, which promote the considered effects. 
Moreover, the long-range magnetic order (due to RKKY 
interaction) is least sensitive to disorder and ME order 
is an intermediate case. Thus we come to the conclusion 
that disorder in the above substances would not adversely 
influence initial ME order. 

It is also seen from the Fig.l that since RKKY inter- 
action is dominating in the problem under consideration, 
the other important factor needed for the ME effect to 
be realized is a sufficient strength of spin-spin (RKKY) 
interaction via free carriers. If T c related to RKKY in- 
teraction is around 150 K, then (see Fig.l) the ME phase 
transition temperature T c me will be of the same magni- 
tude. Moreover, the maximal value of T c me occurs for 
a = 0, i.e. for a two site orientable dipole. For other an- 
gles a and parameters kpd the situation is also not that 
critical, the only "forbidden" values kpd are those where 
PFsd = 0, a = 90° (see Fig.l). This shows that there 
is large space of parameter values where the ME effect 
mediated by free carriers can be realized in the above 
substances. 

To conclude, we have shown that it is quite probable 
to obtain a sufficiently strong ME coupling mediated by 
free charge carriers in mixed valency oxides. Namely, in 
the substances, where RKKY interaction between local 
spins is not small, and the free carrier concentration is 
related to the off-center shift, a substantial ME coupling- 
may be realized. Smaller (but still sensitive) values of 
above coupling occur in the entire domain of kpd values 
except those, where Wsd = (see Fig.l). Our prelimi- 



nary analysis shows, that the effects of disorder influence 
multiferroic behavior in the above substances but do not 
suppress it. The results of quantitative consideration of 
disorder effects leading, for example, to the dependence 
of ME characteristics on the sample history [26| will be 
published elsewhere. 
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